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It is known that classical electromagnetic radiation at a frequency in resonance with 
energy splittings of atoms in a dielectric medium can be described using the classical sine- 
Gordon theory. In this paper we quantize the electromagnetic field and compute some 
quantum effects by using known results from the sine-Gordon quantum field theory. In 
particular, we compute the intensity of spontaneously emitted radition using the thermo- 
dynamic Bethe ansatz with boundary interactions. 
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1. Introduction 



The importance of integrable non-linear partial differential equations in classical non- 
linear optics lias been recognized for some time. The most well-known example concerns 
weakly non-linear dielectric media, where only the first non-linear susceptibility is consid- 
ered important. In this situation, the envelope of the electric field satisfies the non-linear 
Schrodinger equation. This was first understood theoretically by Hasegawa and Tappert|l[]. 
The solitons predicted in |IJ were observed experimentally in [|J. The occurence of these 
classical solitons in common optical fibers promises to revolutionize high-speed telecom- 
munications. 

Our interest in this subject concerns the possibly interesting quantum effects which 
arise when the electromagnetic field is quantized. The latter quantization amounts to 
studying an interacting quantum field theory in the classical variables which satisfy the 
non-linear differential equation. These quantum integrable models have been studied ex- 
tensively, and many of their properties have already been exactly computed. A priori, one 
expects quantum effects to be small for macroscopically large objects such as solitonsB. 
Nevertheless, using known exact results from the quantum non-linear Schrodinger theory, 
quantum effects have been predicted and measured @@] @ li- 
lt is well-known that the non-linear dielectric susceptibilities are enhanced when the 
radiation is in resonance with the energy splitting of quantum states of the atoms of 
the sample. Near resonance, one is no longer in the weakly non-linear regime (higher 
susceptibilities involving higher powers of the electric field become as important) and the 
physics is no longer well described by the non-linear Schrodinger equation. Remarkably, 
as was understood and demonstrated experimentally by McCall and Hahn0, the system 
is well described classically by another famous integrable equation, the sine-Gordon (SG) 
equation. The phenomenon is referred to as 'self-induced transparency'. 

In this paper, we study the quantum effects which arise when electric fields in res- 
onance with a dielectric medium are quantized by using known exact results for the SG 
quantum field theory. In particular, we compute the intensity of spontaneously emitted ra- 
diation by using the thermodynamic Bethe ansatz with boundary interactions^. A different 
but related model is the Dicke Model, which was solved using Bethe ansatz in ||. 

1 In fiber optic systems the soliton consists of a cluster of I0 8 or more photons. 

2 The first few sections of this paper previously appeared in Q . 
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2. Classical Theory 

In this section we review the manner in which the classical sine-Gordon equation arises 
in resonant dielectric media 0. 

We consider electromagnetic radiation of frequency u propagating through a collection 
of atoms, where the frequency ui is in resonance with an energy splitting of the atomic 
states. For simplicity, we suppose each atom is a two state system described by the 
hamiltonian Hq with the following eigenstates: Ho\ipi) = — ^Tiuo\tpi), Ho\ip2) = o^olV^}) 
such that huo is the energy difference of the two states. In the presence of radiation, the 
atomic hamiltonian is 

H atom =H -p-E, (2.1) 

where p = e rj is the electric dipole moment operator. 

We assume the radiation is propagating in the x direction, and E = nE(x,t), where 
n ■ x = 0. The only non-zero matrix elements of the operator p ■ E can be parameterized 
as follows: 

^ 2 \p-E\^)=pE{x,t)e-^, (2.2) 

where p and a are constants which depend on the atom in question. (We have assumed 
spherical symmetry.) It is convenient to introduce the Pauli matrices cr^, and write the 
hamiltonian as 

Hatom = --hLO a 3 - E(x,t)(p 1 a 1 +P2V2), (2.3) 

where p\= p cos a, P2 = P sin a. 

The dynamics of the system is determined by Maxwell's equations, 

{dl - =U 2 ) E{x,t) = J a t 2 P(x,t), (2.4) 

where c 2 = c 2 /eo, eo is the ambient dielectric constant, and P = nP(x,t) is the dipole 
moment per unit volume. The latter polarization can be expressed in terms of the expec- 
tations of the Pauli spin matrices (cr^) = (il)\<Ji\i(j) , where is the atomic wavefunction. 
Namely, 

P = nn(p 1 (a 1 )+p 2 (a 2 )) : (2.5) 
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where n is the number of atoms per unit volume!. Thus, in addition to the Maxwell equa- 
tion ( |2.4[ ), one has dynamical equations for the polarization P(x : t) which are determined 
by Schrodinger's equation for the atom: 



ihd t (<Ti) = (il>\[<Ti,H a t om ]\il>). (2.6) 
The latter can be expressed as 

dt(<?i) =J2^ijkVj (a k ), (2.7) 

where e is the completely antisymmetric tensor with £123 = 1, and 

2E(x,t) 2E(x,t) 
Vl = r Pi, V 2 = ^ p 2 , V3 = u . (2.8) 

To summarize, the dynamics is determined from the coupled equations of motion 
( |2.4|) and (277), wherein the atoms are treated quantum mechanically and the radiation is 



classical. 
Let 

E(x,t) = S(x,t)cos(ujt-kx) (2.9) 

where uj/k = c and S(x, t) is the envelope of the electric field. We will assume the envelope 
is slowly varying in comparison to the harmonic oscillations: dtS <C u)£, d x S kS. In 
this approximation, one finds 



dl-^dA E{x,t)^ 



d x + ^d t )£{x,t) 



sin(u;t-kx). (2.10) 



Let us define (cr»),(a±) as follows 

(an) = (ai) cos(ujt — kx + a) + (a 2) sin(c<;t — kx + a) 

(2.11) 

(aj_) = — (ai) sm(uit — kx + a) + {(T2) cos(a;t — kx + a). 

We make the further approximation that cos 2 (cut — kx + a) terms in the equations of 
motion for dt{(Ji) can be dropped in comparison to cos{ut — kx + a) (and similarly for the 



3 To be more precise, (ai) here represents average over many atoms in a small volume, and is 
thus a continuous field depending on x, t. 



sine terms )□. (This amounts to replacing sin 2 (ut — kx + a) : cos 2 (ut — kx + a) by 1/2) . One 
then finds 

d t (cr\\) = (u-uo) (a±) (2.12a) 

d t (a ± ) = -^lp (a 3 ) + (wo - W )<CT||> (2.126) 

W = ^W- (2.12c) 

Finally, equations ( p.4f) and ( |2.1C| ), upon making an approximation analagous to the slowly 
varying envelope on the RHS of (|2.4j), lead to 

(d x + ld t )£(x,t) 
c 



2tt 

sin(u;t — kx) = —cnivp (sin (cut — kx)(a±) — cos(ivt — kx)(a^)j 



(2.13) 

In order to solve these equations, note that d t (J2i( a i) ( a i)) = 0- Thus, if the atoms 
start out in their ground state l^i), with (0-3) = 1, one has J2i( a i) 2 = 1 f° r an time. 
On resonance, when uj = ujq, eq. ( |2.12q| ) implies (cry) = for all time. The constraint 



Y2i( a i) 2 = 1 can be imposed with the following parameterization: 

(ffi) = sin(/3 cl <j)(x, t)), (a 3 ) = cos(P cl (P(x, t)). (2.14) 

The parameter (3 c i is arbitrary at this stage, but will be fixed in the next section. Equations 
( p. 126, cD now imply 



d t <t> = --^-£(x,t). (2.15) 
Inserting ( |2.15| ) into ( |2.13|) and defining 

x' = 2x-ct, (2.16) 

one obtains the SG equation: 

(a 2 -c 2 a 2 ,)^ = -£sin(^), (2.17) 



where 

u 2 = - 

he 



^ = 2 -^. (2.18) 



It can be shown in perturbation theory that this is a good approximation at or near resonance. 
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3. The Quantum Spectrum 

In this section we proceed to quantize the electromagnetic field. In order to do this, 
the SG field <p must be properly normalized such that the energy of soliton solutions 
corresponds to the physical energy; this amounts to properly fixing the constant /3 c i. 

The action which gives the classical SG equation of motion is 

S so = \f dx'dt QW) 2 - - (d x ,cp) 2 + ^ co S ((3 cl( p)^ . (3.1) 
On the other hand, the properly normalized Maxwell action is 

^Maxwell = ^ J d s xdt (-If^F^ + ...J = ^ J d 3 xdt \{d t A) 2 + .... (3.2) 

where the vector potential is A = nA and as usual E = —^dtA. (One can show that 
Aq can be set to zero.) One should be able to show directly that the complete Maxwell 
action, upon making the appropriate approximations of the last section, reduces to the 
sine-Gordon action, since we have already shown this to be the case at the level of the 
equations of motion. Assuming this to be the case, in order to properly normalize the field 
0, one needs only compare the kinetic terms in (|3.1|) and ( |3.2| ). The dimensional reduction 
is made by assuming simply that A is independent of y, z and f dydz = A, where A is 
an effective cross-sectional area perpendicular to the direction of propagation (as in the 
cross-sectional area of a fiber) . From ( |2.15| ) one has 



'Maxwell 



A {^y^j J dxdt ^(d t (p) 2 cos 2 (ut -kx) + .. 



(3.3) 



In the above equation, we make an approximation analagous to the one made in arriving 
at (|2.12|) and replace cos 2 (^t — kx) by 1/2. Noting that dxdt = dx'dt /2, upon comparing 
with Q, one fixes (3 2 cl = Ap 2 /(Ah 2 c). 

Finally, as is conventionally done in the quantum SG literature, we rescale <fi — ► Vh<fi, 
so that SsG/h takes the form (|3.1|) with (3 c i replaced by /3, where 

The constant j3 is the conventional dimensionless coupling constant in the quantum SG 
theory, and is allowed to be in the range < /3 2 < 8n, where (3 2 = 8tc corresponds to 
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a phase transition [ TO] . The limit where the radiation is classical but the atom is still 
quantum mechanical corresponds to the limit f3 — > 0. 

The classical soliton solutions to the SG equation are characterized by a topological 
charge T = ±1, where 

T = ^-(<P( X ' = oo) - (j){x = -oo)). 

Solitons of either charge correspond to solutions where at fixed x the atoms in the far 
past are in their ground state, and are all in their excited state at some intermediate time: 
( (T 3)t=±oo = 1- These classical solitons have been observed experimentally in 0. What 
distinguishes solitons (T = 1) from antisolitons (T = — 1) is the sign of the envelope of the 
electric field. From the known classical soliton solutions and ( |2.15| ) one finds for T = ±1, 



^^vKHv^J) ■ (3 - 5) 

Thus the electric fields for the soliton versus the antisoliton are out of phase by n. 

The particle spectrum of the quantum field theory is knowni to consist of soliton and 
anti-soliton of mass m s , and a spectrum of breathers of mass 

TV) StT 

m n = 2m s sin— , n = l,2, .... < — , (3.6) 
16 7 

where 

7 = -f^- (3-7) 

The n-th breather disappears from the spectrum when 87r/7 = n; at this threshold the 
mass of this breather is twice the soliton mass. The n-th breather can thus be considered 
as a bound state of n soliton/antisoliton pairs. The n = 1 breather is known to correspond 
to the particle associated with the sine-Gordon field </> itself. 

From ( |3.2|) ( p73D , one sees that the sine-Gordon scalar field is essentially the photon 
vector potential, thus the n = 1 breather may be thought of as a 'envelope photon', which 
is massive due to the dielectric properties of the medium. Of course, this photon has little 
to do with the real free photons of energy Tiuj. The n-th breather state can then be viewed 
as an n-envelope-photon bound state. Under this heuristic assignment of photon number 
to the spectrum, the soliton has photon number 1/2. In specific physical situations, f3 is 



5 The literature on the quantum SG theory is extensive. See for example [11], and references 
therein. 
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very small (see below), thus mi << m s , which indicates the strong binding energy of a 
pair of solitons to form an envelope photon. 

The soliton mass is directly related to the parameter fi in ( |2.18j ). In the quantum 
theory, short distance singularities are removed by suitably normal ordering the cos (/?</>) 
potential |10||: 

-pr cos(/30) — > A : cos(/30) : 
The anomalous scaling dimension of the operator : cos(/30) : is /3 2 /47r, so that A has mass 



dimension 2(1 — |^). Therefore 



m s oc (p) 



V(i-fc) 



(3.8) 



This allows us to obtain quantum corrections to the frequency and density dependence of 
the soliton mass. Since \i oc y/um, one finds 



P 2 



m s oc y/um I 1 + ^— log(um) + 0(f3 4 ] 
Ibn 



The exact relation between m s and \x is known W\: 

m s = c (/3)m 7//32 A- 7/8w , 
where A is an ultraviolet cutoff, and 



m = 



11% 

c z 



2r( 7 /167r) 



n 



r(i-fe: 



2/3 2 r( 7 2 /87r/3 2 ) 



T/2/3 2 



As /3 — > 0, one obtains the well known classical expression, 

(/3-0). 



8m 



(3.9) 



(3.10) 



(3.11) 



(3.12) 



From ( |3.6| ) , one sees that for small (3, m is approximately the mass of the lowest breather. 

The magnitude of the quantum corrections to classical results is determined by the 
parameter /3 2 /87r. In the original experiment of McCall and Hahn, a ruby sample was 
used withp = 4.8 x 10~ 21 in cgs units, which corresponds to (3 2 = 10~ 23 y/eo/A. The ruby 
rod had a cross-sectional area of A ~ 1cm 2 . Taking the density n to be that of the Cr +S 
doping atoms in the ruby, then n w 10 19 /cm 3 . For ui w 10 15 s _1 , one finds m w lO -37 ^, 
and m s ~ 10 _13 5r. This large hierarchy of mass scales between the soliton and the lowest 
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breather is due to the smallness of f3. Thus the experiment performed by McCall and 
Hahn is in the extreme quasi-classical limit, with an extensive spectrum of light and heavy 
breathers and relatively heavy solitons with macroscopic masses. In the quantum theory, 
the classical soliton observed in for example 0, becomes a fundamental particle, i.e. it is 
not to be thought of as a coherent state of photons, due to the fact that the soliton carries 
non-zero conserved topological charge. 

Since in the limit /3 — ► 0, m s = 8m/ P 2 , from one can express /3 2 /87r in terms of 



the classical soliton mass: 



Thus, 



1 ' A — (3.14) 



871" iV 7 V-t^atom 



where iV 7 = %uj/m s c 2 roughly corresponds classically to the number of photons that 
comprise the soliton, A c is the Compton wavelength of the soliton, and L atom = 1/nA 
is the inter-atomic spacing. The above equation summarizes where one expects quantum 
effects to be important: when the soliton is comprised of small numbers of photons, or 
when the Compton wavelength is large compared to the space between the atoms. 

The classical scattering matrix for the solitons has been computed by comparing N- 
soliton solutions in the far past and far future lfL3f . The exact quantum S- matrix is also 
known |TlJ| . Since ft 2 /8n is small here, the quantum corrections to classical scattering are 
most easily determined by incorporating the one-loop corrections to the classical scattering, 
which amounts to the replacement f3 2 — > 7. The necessary formulas can be found in the 
above papers. 

The most interesting aspect of the quantum scattering of solitons is that the so- 
called reflection amplitude is a purely quantum effect, analagous to barrier penetration 



rjj] . Namely, one considers an in-state consisting of a soliton of momentum pi and an 
antisoliton of momentum P2 which scatters into an out state where the momenta p± and P2 
are interchanged. The T = +1 soliton is thus reflected back with momentum equal to that 
of the incoming T = — 1 antisoliton. In the semi-classical approximation, this reflection 
amplitude is 

S R (9) = \ [e^h _ {j e -^\o\ s{6) (3 . 15) 



where 



exp - 

,7 Jo 



drj log 



e & - %r > + 1 
e 9 + e~ ir > 



(3.16) 



and 9 = 9i — 62, where pi 5 2 = M s sinh^i^- Since I/7 ~ 1//3 2 is very large, the oscillatory 
factor exp(167r 2 z/7) — 1 in (|3.15| ) will make the detection of reflection processes difficult, 
since even in a small range of /3 2 , this factor averages to zero. 



4. Spontaneous Emission 

In this section we will consider the spontaneous emission of radiation from a collection 
of atoms in their excited state. We suppose that all the atoms have been 'pumped' by some 
means to their excited state at t = 0. In the coordinate system (x',t'), defined by ( |2.16j ) 
and t' = t, the atoms are along the x' axis at t = 0. If is the initial state corresponding 
to all atoms being in their excited state, then 

(#r 3 |V> = -l = oos/ty, (t = 0). (4.1) 

Thus, the initial state corresponds to the initial constant field configuration 

<P = <Po = ^- (t = 0). (4.2) 

All atoms initially in the ground state corresponds to </>o = 0. The subsequent time 
development of the system thus amounts to an initial value problem in the SG quantum 
field theory. 

The initial condition ( [4.2|) may be imposed by adding to the SG action the boundary 
term: 

^boundary = 9 J dx COS {(3{(j) ~ (f) )/2) . (4.3) 

In the limit g — > oo, this enforces the boundary condition <p = <pQ. It is known that the 
boundary interaction ( |4.3| ) preserves the integr ability of the SG theory [|T5||. Ghoshal and 
Zamolodchikov have developed a comprehensive approach to studying integrable quantum 
field theory with boundary interactions, and in particular have solved the boundary SG 
theory. The limit g — > oo is referred to as the fixed boundary condition in [j!5| , whereas 
g = is called the free boundary condition. As shown there, information about the 
boundary condition is concentrated in a boundary state \B) at t = 0, with the general 
form 

\B) = exp d6K ab (6)A a (-6)A b (6)^) |0> (4.4) 
9 



where A a are creation operators of particles of type a with rapidity 9, and 



K ab (9) = R^/2-9), (4.5) 

where R{9) are the S-matrices for reflection off the boundary, and a is the charge conjugate 
of a. The rapidity parameterizes energy and momentum as follows, 

E' = m cosh 9, P' = m sinh 9, (4-6) 

where E',P' are energy and momentum in the coordinate system (x',f). (We have set 
h = c = 1.) In our application, since the boundary is a spacial axis at t = 0, K ab (9) 
represents an amplitude for emission of particle pairs at the boundary. 

In order to compute the intensity, one should first study a partition function with the 
relevant boundary conditions. In euclidean space, one may consider the partition functions 

Z(R) = (B'\e- HR \B), (4.7) 

where H is the hamiltonian, and (B'\ is a conjugate boundary state. Taking \B) to be 
the boundary state for the fixed boundary condition with </>o given in ( |4.2| ), and (B'\ to 
be the boundary state for the free boundary condition provides information on the time 
evolution of a state where the atoms are all initially excited, and after a time interval R 
the final state has no constraints on <fi. 

To simplify the initial discussion, let us assume there is only one type of particle 
in the spectrum, with bulk two-particle S-matrix S(9), and boundary S-matrices K(9) 
and K'(9) for the boundary states \B) and \B') respectively. The function Z(R) has the 
interpretation as a partition function on a strip of width R with boundary conditions \B) 
at both ends. As shown in |1(J, if one uses the expression (|4.4| ) for \B), then one encounters 



divergences in Z(R) due to the fact that the formula (|4.4j ) is only valid in infinite volume. 
We thus regulate these divergences by taking the boundary at t = to have a large but 
finite length L. (L corresponds to the length of the sample.) Then Z(R) is characterized 
by the thermodynamic Bethe ansatz (TBA)[^6|. These TBA equations generalize the bulk 
equations with no boundaries derived in jTTj . One introduces a density of pairs of particles 
per unit length p{9) and similarly a density of holes p h {9). Due to the finite volume L, the 
momenta are quantized in a way that involves the bulk S-matrix: 

p{9) + p h {9) = B coshfl + i (<& * p) (9), (4.8) 

ZTT ZTT 
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where 

$(d) = -id e logS(d), (4.9) 

and * denotes convolution, 

/oo 
dd'${6-6')p{d'). (4.10) 
-oo 

Defining 



- = e £ , (4.11) 
P 



then e satisfies the integral equation 



e{6) = 2mRcosh6 - log(K(K')*) - * log (l + e" £ ) . (4.12) 

Z7T 

The partition function is then given by the formula 

log Z(R) = / d6> cosh 6> log ( 1 + e- eW ) . (4.13) 

4tt J_ 00 V / 

One can now compute the intensity, i.e. the energy flux per unit time of the emitted 
radiation. This radiation will consist of waves traveling to the right and left. Above, the 
SG theory describes the envelope for the right moving wave ( |2.9|) , therefore we compute 
the intensity in the right-moving radiation. The TBA analysis in |16| assumes periodic 
boundary conditions in the spacial direction (cylindrical geometry), however we consider 
the intensity computed below to be a good model of the measureable intensity at one end 
of the sample. 

In the slowly varying envelope approximation, the Pointing vector is 



S = ^^£ 2 (x, t) cos 2 (ujt - kx)x. (4.14) 

47T 

Averaging over periods of the coswt oscillations, one has that 

/ = ^e„v, (4.15) 

where / en v is the intensity of the envelope. J env can be computed from the particle picture 
of the above TBA. Let (E,P,v), (E',P',v') be the energy, momentum and velocity of a 
single particle in the (x, t) and (x ; , t') coordinate systems respectively. One has 

P = 2P', E = E'-cP', v = (v' + c)/2. (4.16) 



11 



If E', P' are given as in (|4.6| ), then E = me" , and v = (1 + tanh#)/2. 

In the coordinate system (x', £'), pairs of particles are emitted with rapidity 9 and —9. 
For both particles in the pair, v(9) > 0, thus 



/•OO 

/env = / d9p(9) (E(9)v(9) + E(-9)v(-9)) 
Jo 



, .x , (4-17) 



1 d9p(9) 



cosh 6* 



(We have reinstated c, h; p has units of (length) -1 and satisfies the equation (|4.8|) with 
m — > rric/%. ) 

In the SG case, for general /?, the above TBA equations are very complicated due to 
the non-diagonal pieces (reflection pieces) of the soliton S-matrix. To simplify the situation, 
we recall that (3 2 /8tv is very small in physical applications, and is thus well approximated 
by 

3 2 1 

SF = wTT (418) 

for some integer N . Namely, iV is taken as the integer which best approximates (|4.18|) . At 



the points ( f4.18| ), the reflection amplitudes vanish and the scattering becomes diagonal. As 



explained above, for the problem of spontaneous emission, one should take the boundary 
state (B'\ to be for the free boundary condition. This leads to the problem that the solitons, 
though they scatter diagonally off the fixed boundary condition, scatter off-diagonally from 
the free boundary condition, and it is not known how to formulate TBA equations in this 
situation. Instead we solve a different problem corresponding to an initial state wherein 
all the atoms are in their excited state, and a final state where all the atoms are in their 
ground state. Both of these are fixed boundary conditions with diagonal soliton scattering 
off the boundary. 

The TBA equations become the following. Introduce densities of particles and holes 
per unit length for each type of particle p a (9), p } ^{9) where a runs over breather and 
soliton indices: a G {1, 2, N — 1, +, — }, where ± refers to the soliton and anti-soliton. 
Introducing 

^=e £ % (4.19) 

Pa 



one has 



Pa + p h a = T ^ COSh 9 + ^ ab * Pb (4.20) 
12 



v « ~ T~ S ® ab * log ^ + e " £b ) ( 421 ) 



2tt 
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L /"°° 

logZ(i2) = — > / d# m a cosh#log(l + e _£a ) (4.22) 

47T / 

, env = E ^r^ W _i_, (4.23) 



where 

u a = 2m a Rcosh9- log{K a {K' a )*), (4.24) 

$ ab = -id 9 logS ab (e), (4.25) 

where is the 2-body S-matrix for particles of type a, b and K a are the boundary S- 
matrices for fixed boundary condition 0o =7r /(3, and are those for the fixed boundary 
condition <po = 0. 

The SG scattering theory at these reflectionless points is equivalent to the minimal 
scattering theory associated to the root system of the Dn+i affine Lie algebra JT5| ]T5 



There is a remarkable universal form of the above TB A equations [ 20 ] [^T 



S a = "a - ^ Yl Nab ® * ^ b ~ L °S i 1 + e£b )) > ( 4 - 26 ) 

where $ is the universal kernel: 

* <9) = SOS- (42T) 

and iV^ is the incidence matrix for the Dynkin diagram below. Using similar manipulations 
as in |20|j ||21| , one can obtain the universal form of the integral equation for p a , p^: 

Z7T z — ' 

6 



where 



777 

^cosh#. (4.29) 

2/T 



The boundary S-matrices for the fixed boundary conditions can be found in [ITS 



Specializing to the reflectionless points, and using some Gamma function identities, one 
finds the following!. We give below the R reflection amplitudes; the K functions that 



6 The denominator of (5.23) in |HJ should be n(z,7r/2)n(-z,7r/2)n(a:, -7r/2)II(-x, — 7r/2), 
and (5.24) should read a(x,u)a(x, —u) = cos 2 x [cos(x + Xu) cos(x — Ait)] -1 . 
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N-2 



Figure 1. Dynkin diagram for D^+i- 




appear in the integral equations above are given by K a (9) = R a (iir /2 — 6). For the 
solitons, one has 

JV-l 



R ± (u) = (-1) 



where u = —iO, £ = -S-(f>o, and 



N- 



. lC os(e±iV«)^ , A J fr isin (f + m) 



B 

N-l 



= J] 



COS 



2 AT 



COs£ 



7r(l+2fc) 
4N 



n - 



LA sinf- - —) ' 



(4.30) 



cos 



27V 



7r(l+2fc) 
AN 



For the n-th breather one has 



=o cos ( ^ + ^ (1 + 2fc) - li/2 ) cos ( ^ - 



7r(l+2fc) 
47V 



(4.31) 



R (n) (u) = , 1)n+1 cos (f + j|) cos (| - f - f|) sin (| + I ) 
l " "' i!)'-(f'T'f)^(i T) 

sin( W +ff)cos 2 (f -f 



(*), 



cos(f 

n-l 

X H^-S)cos 2 (f + f + &) 



(4.32) 



(4.33) 



4 2n V)=n 



'sinu - cos (|p - (/ - l/2)f ) ^ ^sin u - cos (£ + (/- l/2)f ) 

^sinw + cos (I + (/- l/2)f ) 



y sin u + cos 1/2) 
n = 1, 2, • • • < JV/2 



jR (2n-i) , x = cos(£/iV) - sin-u yr 

cos(£/iV) + sin u 11 sin „ + cos (| _ £ j 

Ar + 1 

n= 1,2,...,< — — . 



7 7 



sin it — cos ( — I sinw — cos ( + 



sin u + cos 



\N ^ N J 



(4.34) 
(4.35) 
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5. Concluding Remarks 

We have shown how quantized radiation interacting with two-level atoms leads to the 
quantum sine-Gordon theory in the deep semi-classical regime. As expected, this implies 
quantum corrections are quite small. We have not considered here the kinds of quantum 
effects which were studied previously for the non-linear Schrodinger case, which involve 
taking into account coherent state initial conditions of the radiation in real systems and 
the squeezing of such states. For the resonant situation considered here, the topological 
charge of the sine-Gordon solitons presents difficulties for constructing analagous coherent 
states of photons. 

An interesting theoretical question which arises in this work concerns the classical 
limit, ft — > 0, of the thermodynamic quantities in the SG theory. As suggested above, 
one can approach this limit by hopping along the refectionless points (|4.18|). Then the 



problem becomes that of the N goes to infinity limit of the -Dat+i TBA equations. Since 
this corresponds to a certain classical limit of the SG theory, this leads us to believe this 
limit is meaningful. The interesting limit from a physical point of view is to take ft — > 0, 
keeping m s , which tends to 8////3^jC 2 , fixed. This is reminiscent of the large iV studies 
of the ground state energy for the SU(N) sigma models (principal chiral model), where 
analytic expressions were obtained in terms of Bessel functions p3|. Another interesting 
question concerns how sensitive the classical limit is to the way in which ft tends to zero. 
I.e. is the limit different if one does not hop along the reflectionless points? We hope to 
address these issues in a future publication. 
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